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Recent experiments on the amorphous magnetic semiconductor Gd^Sii-a,, Phys. Rev. Lett. 77, 
4652 (1996), ibid 83, 2266 (1999), ibid 84, 5411 (2000), ibid 85, 848 (2000), have revealed an 
insulator-metal transition (i-m-t), as a function of doping and magnetic field, a spin glass state 
at low temperature, and colossal magnetoresistance close to the i-m-t. There are also signatures of 
strong electron-electron interaction close to the i-m-t. Motivated by these results we examine the role 
of doped magnetic moments in a strongly disordered electron system. In this paper we study a model 
of electrons coupled to structural disorder and (classical) magnetic moments, through an essentially 
exact combination of spin Monte Carlo and fermion exact diagonalisation. Our preliminary results, 
ignoring electron-electron interactions, highlights the interplay of structural and magnetic 'disorder' 
which is primarily responsible for the observed features in magnetism and transport. 



I. INTRODUCTION 

A. Insulator-metal transition in doped 
semiconductors 

The role of disorder in electronic systems is an endur- 
ing theme in condensed matter physicsEm. Structural 
disorder leads to localisation of electronic states and, in 
a three dimensional system, all electronic states would 
be localised^ if the disorder, 'A' say, were greater than a 
critical value (A c ). For weaker disorder, A < A c , only 
states beyond an energy e c (A) of the band center, i.e in 
the band tails, are localised. A system is metallic or 
insulating depending on whether the Fermi level, ep , is 
above or below this 'mobility edge' e c . Variation in elec- 
tron density, and hence ep, can drive a system through 
the i-m-t. This is the basic scenario for the 'Anderson 
transition' in a non interacting electron systemtfu. 

Doped semiconductors^] have served as a laboratory for 
studying disorder effects since they allow systematic con- 
trol of the carrier density. The i-m-t has been studied in 



most famous example be- 
(Si:P), and less extensively 
. In experimental systems, 



doped crystalline systems, tl 
ing phosphorus doped silicon! 
in amorphous semiconductor; 

however, the phenomena of localisation is complicated by 
k> \ electron-electron (e-e) interaction .-effects which become 
5_J ■ very pronounced close to the i-m-ta. There are clear sig- 
C3 ' natures of interaction effects in the density of states and 
the conductivity in amorphous systems close to the i-m-t. 
The interplay of disorder and interaction effects near the 
i-m-t is not completely understoodO'tTj and continues to 
be an area of active research. 

The introduction of magnetic moments in a 'disor- 
dered' system brings in new phenomena. The affects are 
well understood in crystalline semiconductorsEJI13, where 
the disorder is weak and the relevant regime is of low car- 
rier density. In these systems, electrons are trapped in 



the potential fluctuations and polarise the magnetic mo- 
ments in their neighbourhood. This leads to a gain in 
exchange energy, and tends to enhance the localisation 
of the carrier. The net localising effect in these systems, 
therefore, arises from a combination of (£) on site disor- 
der; a 'single particle' effect, and (ii) 'bound state' for- 
mation between the electron and the magnetic moments; 
a 'many body' effect. The composite object, a (trapped) 
electron and the polarised spins in its vicinity, is called 
a 'spin polaron'. Since the i-m-t in these systems occurs 
at a carrier density ~ 10 -5 (per unit cell), the picture 
of non overlapping, uncorrelated, spin polarons is ade- 
quate. The interesting difference between these systems 
and their non magnetic counterparts lies in the response 
to a magnetic field. An applied field globally aligns the 
magnetic moments, reduces the polaron binding energy, 
and tunes the 'magnetic component' of localisation. The 
change in activation energy (or mobility) aLihe carri- 
ers can lead to enormous magnetoresistance!2't3. These 
systems have been known to exhibit 'colossal magnetore- 
sistance' (cmr) much before such effects were observed in 
the manganites. 



Recent experimentsE-3~li3 on amorphous (a)-GdSi re- 
veal that doping magnetic moments in an amorphous 
system combines the richness, and complication, of the 
traditional i-m-t in amorphous semiconductors^ with the 
physics of 'cmr'. In addition a disordered magnetic state, 
a spin glass, with rather unusual properties emerges at 
low temperature. This combination of an i-m-t, in a 
high carrier density system, with cmr, and a spin glass 
ground state is probably unique to GdSi. We will dis- 
cuss these experiments in the next section. The essential 
phenomena in these systems seems to be: (i) electron 
localisation due to the strong disorder in the amorphous 
structure, (ii) indirect exchange interaction between the 
doped moments, mediated by the electrons, leading to 
a spin glass state, (Hi) feedback of the spin background 
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on the electronic system, via 'spin disorder', tending to 
enhance electron localisation, (iv) strong e-e interactions 
in the vicinity of the i-m-t, showing up as a correlation 
gap in the density of states, a y/T dependence of the low 
temperature conductivity, and 'local moment formation' 
in the electron system, and, (v) control of the 'spin dis- 
order' by a magnetic field, h, leading to large negative 
magnetoresistance when the electron density is close to 
the critical density. 

Unlike the case of crystalline magnetic semiconduc- 
tors, wheats— .the picture of isolated 'bound magnetic 
polarons'tStS seems to suffice, we do not have an un- 
derstanding of a high density of electrons in the back- 
ground of strong disorder, interacting with randomly lo- 
cated magnetic moments and with each other. Scenarios 
in terms of an Anderson transition, or isolated polarons, 
or model spin glasses, do not suffice, and the interplay of 
these effects is what we study in this paper. 

While we are motivated by the phenomena in GdSi, 
which we take to be the prototype amorphous magnetic 
semiconductor, our focus in this paper differs from the 
experiments in two respects, (i) In the experiments, dis- 
cussed later, the tunable parameters are the electron den- 
sity (n e i) and the density of doped magnetic moments 
(n sp ). Carrier density and 'moment density', however, 
are not independent variables since both come from the 
doped Gd atoms. Thus, n sp ~ n e i and most of the data 
is for n sp = x ~ x c , with variations of a few percent. 
The 'disorder' (A) due to the amorphous structure, or 
the coupling (J 1 ) between the electron and doped mo- 
ments cannot be varied experimentally. We will discuss 
a model which is expected to describe GdSi (see later) 
but not restrict ourself to parameters specific to the real 
system right away. The model has a large parameter 
space and we will explore it gradually to locate some of 
the phenomena, i-m-t, spin glass and 'cmr', mentioned 
earlier. We will present our results specific to GdSi in a 
separate publicationta. (ii) The effects in GdSi arise from 
a combination of structural disorder, electron-spin inter- 
action, and e-e interaction. All of them are important, 
in varying degrees, for the phenomena observed. We will 
study a model problem of electrons coupled to structural 
disorder and magnetic moments, and ignore e-e interac- 
tions in the present discussion. They are important, as 
we will see in a detailed review of the experiments, but 
we have a non trivial problem even without such interac- 
tions. As far as we know even this simplified model has 
not been explored. We will estimate and quantify-the e-e 
interaction effects in our results specific to GdSiia. 

Here is the outline of the text. In the next section 
we discuss the experimental results on amorphous (a)- 
GdSi in some detail, to highlight the effects of doping 
magnetic moments into a disordered system. Following 
that we define a model which we believe contains most 
of the relevant physics. We then discuss the approxima- 
tions involved in 'solving' for the properties of this model 
and present results on the magnetic properties, thermo- 
dynamics, and some simple limits for the conductivity. 



We conclude by indicating how electron-electron interac- 
tion effects can be included, approximately, within our 
scheme. This would recover some of the features in the 
experimental data which are not accessible in our 'non- 
interacting' model. 



B. Magnetic moments in an amorphous background: 
experiments 

Treating GdSi as the 'model' amorphous magnetic 
semiconductor we review some of the effects which have 
been experimentally observed. The measurements have 
been madeEjilfl on a-GdSi and simultaneously on the 
non-magnetic analog a-YSi, to clarify the role of mag- 
netic moments on the i-m-t. The results are broadly 
on (i) the conductivity, er, and the magnetoresistance, 
(ii) spectral properties/density of states (DOS), probed 
through tunneling, (Hi) thermodynamic properties: spe- 
cific heat, CV, and entropy, S, and (iv) magnetic prop- 
erties: the linear response susceptibility x(T) and the 
magnetisation m(h, T). ._. 

Let us begin by with the conductivityc3. (a). Both 
Y^Sii-a, and Gd^Sii-a, show an insulator-metal transi- 
tion as the doping, x, is increased across a critical value 
x c . The critical doping required (x c ~ 14%) is slightly 
greater in GdSi compared to YSi. (b). For both GdSi and 
YSi, in the doping range 0.17 > x > 0.11, the conduc- 
tivity has da/dT > 0, i.e dp/dT < 0, at all T. This is 
true even of the systems which are 'metallic' at low tem- 
perature, (c). For the metallic samples the low temper- 
ature conductivity can be fitted to o~(T) ~ Co + o~\\JT. 
(To and <ti are constants and the yT term arises from 
Coulomb effects in a disordered systemEl. On the insu- 
lating side o~(T) ~ e ~v /T °/ T ; indicating variable range 
hopping (VRH) in the presence of a soft Coulomb gapEEl. 

(d) . All Gd based systems, in the doping range studied, 
show pronounced negative magnetoresistance (MR). The 
MR increases on reducing x (towards x c ), on decreas- 
ing T, and on increasing h. YSi samples of comparable 
composition show a significantly smaller, positive, MR. 

(e) . The field dependence of conductivity is seen to arise 
principally from variations in o-q in the metallic phasecll, 
and from the variation of Tq in the insulatorn3. So, the 
dominant temperature dependence of the conductivity in 
the metallic phase arises from Coulomb effects, driving 
the \ff, but the MR arises almost entirely from varia- 
tions in the T = conductivity. (/). For GdSi samples 
with x < x c an applied magnetic field can drive the sys- 
tem metallic. Near this field driven i-m-t, er varies as 
(h — h c ), h c being the field at which the i-m-t occurs. 
The effects, (a) and (d) — (/), hint at the crucial role of 
magnetic moments in determining electronic transport. 

The low energy DOS has been studiedlia, from mea- 
surement of the tunneling conductance, across the field 
driven i-m-t. These measurements provide direct in- 
formation on the evolution of the correlation gapB (dip 
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in DOS at ep) riear the metal-insulator transition and 
the Coulomb gapEJ in the insulating phase. The princi- 
pal conclusion from these measurements is that, for the 
field tuned transition, the density of states at the Fermi 
level, N(0), varies as (h — h c ) 2 while ao(h) ~ (h — h c ). 
Therefore, the critical behaviour near the i-m-t follows 
N(0) oc a 2 , with both the DOS and the d.c conductivity 
vanishing at the transition. 

The difference in the thermodynamic properties of 
magnetic and non magnetic systems again emphasises 
the role of the doped moments. Even for disordered 'non- 
magnetic' systems close to the metal-insulator transition 
the specific heat has a contribution arising from 'local 
moments 'E3. This arises because electrons can behave 
like localised S — 1/2 objec±s.in the presence of disorder 
and strong e-e interactionsEil. However, after 'subtract- 
ing out' the specific heat of the non magnetic analog, one 
would expect the 'high temperature' entropy of the mag- 
netic system to correspond to S°(oo) ~ n sp log(2S + 1). 
Results on GdSi differ from this: (a). Specific heat mea- 
surements indicate that at 'high' temperature S mag (T) = 

J T dT'Cy a9 /T' is larger!! than S°(oo) by approximately 
50%! The total magnetic entropy tends to saturate, at 
this (larger) value by T ~ 60 — 70 K. A simple minded 
estimated suggests that the additional entropy can be 
accounted for by ~ 2 spin 1/2 moments (conduction elec- 
trons) for each Gd moment, (b). Cy aa should vanish as 
T — > 0, however, down to 5 K, it shows no sign of a 
downturn^. This suggests that the peak in Cy ag occurs 
somewhere below 5 K. At high temperature the magnetic 
and non magnetic systems have the same Cy. 

Finally, the magnetism. The doped Gd atoms possess 
a moment J = S = 7/2, arising from the half- filled / 
shell. The coupling between the randomly located mo- 
ments arises primarily through mediation by the conduc- 
tion electrons. The overall magnetic effects also have 
contributions from the 'local moments' in the conduction 
electron system, alluded to in the previous paragraph. 
The basic observations on the magnetism are: (a). Mea- 
surements at high fields, h ~ 1 T, show a magnetisation 
growing as ~ 1/T, i.e, free moment like response. How- 
ever, measurements of x at low field (h ~ 10 mT), re- 
veals a transition!]] to a spin glass state at a temperature 
Tf ~ 1 — 6 K depending on doping, (b). The freezing tem- 
perature increases from Tf ~ 1 K at x = 0.04 to T t <> i 6 
K at x ~ 0.20. The observation of the characteristicE3~E3 
(logarithmic) frequency dependent shift of Tf confirms 
that these are indeed spin glasses, (c). There is a distinct 
difference between the field cooled (FC) and zero field 
cooled (ZFC) susceptibilities x(T); the FC susceptibility 
saturates to a higher value as T — > 0. (d). Fitting a Curie 
form to the susceptibility, x( x > T) = A(x)/(T ~ 9(x)), re- 
veals that the effective moment, /i e //, which can be ex- 
tracted from A, varies significantly with x. Naively, this 
should have been just yj S(S + 1), independent of x. It 
is close to this expected (free moment) value for x m x c 
and falls off on both the metallic and insulating sides. 



(e). The spin glass freezing does not have any signature 
in the temperature dependence of transport properties. 
The freezing itself is eliminated by fields h > 0.1T. 

These results suggest that the magnetic state, which 
crucially affects electronic transport, is itself determined 
by electron spin interaction, the background disorder, 
and e-e effects. 

Let us abstract our lessons from these effects, to start 
on a theory. The primary effect in amorphous magnetic 
systems is electron localisation due to disorder, either 
due to randomness in the amorphous structure or dis- 
order in the spin configuration. To a first approxima- 
tion, the i-m-t is driven by varying the mobility edge 
e c across ep. However, as the data indicates, the tran- 
sition is not a simple 'Anderson transition', driven by 
structural and spin disorder, in a non-interacting electron 
system. Electron-electron interaction effects are clearly 
in evidence in (i) the temperature dependence of c(T), 
(ii) low energy features in the DOS, N(e), and, possibly 
also in (Hi) the 'excess entropy', at high temperature, 
and (iv) the non monotonic variation in fJ. e ff across the 
i-m-t. 

The transition, however, is not driven by interaction ef- 
fects per se. Interactions become relevant near the transi- 
tion, there the effects of disorder and interactions cannot 
be deconvolved, but most of the phenomena occur be- 
cause the system is strongly disordered. Since we will 
study a model without e-e effects, note that many of the 
experimental results can be understood, at least qualita- 
tively, even without such interaction. These include (i) 
the existence of an i-m-t, from variation in the 'effective 
disorder' seen by the electron, (ii) a magnetically disor- 
dered ground state and the spin glass signatures, and, 
(in) the large MR from the field tuning of spin disorder. 

This problem, even without Coulomb effects, is non 
trivial because: (a) there is no simple method for access- 
ing transport properties in the regime we are interested 
in (even if the 'disorder' were completely specified), and 
(b) the 'spin disorder' is not specified, the magnetic state 
itself depends intimately on the electronic spectrum and 
wavef unctions. This coupling is at the heart of the prob- 
lem. 



II. THE MODEL 

The model for amorphous magnetic semiconductors 
would have the form 

H = '^2(Uj, C rcl a Cj cr + h.c) - ^ /m, + j'y^ cr v .S u 

ij i v 

■ X! 1 </"<"< C 1 ) 

ij 

where the label i,j etc refer to the (non periodic) atomic 
locations Ri,Rj etc. The labels v refer to some set of 
positions {iv}, say, where the magnetic ions are located. 
The ta would have a distribution since the structure is 
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non crystalline. The Vij stand for the Coulomb inter- 
action, ov is the electron spin operator and we have 
assumed the local electron-spin coupling, J', to be site 
independent. 

A simpler variant of this model, presumably with sim- 
ilar physics, is the following: 

H = -t 2J (cJo-'V + h.c) + y^(ej - p)n lcr 

+ J 'Y1 <J ' J - Sv + ZZ V i] n * n 3 ( 2 ) 
v ij 

Here the i, j refer to sites on a periodic structure and the 
t now refers to nearest neighbour hopping on that lattice. 
The 'amorphous' nature is incorporated via the random 
on site potential, e,, and the sites {iv} are some subset 
of the lattice points {R;}. The parameters of the theory 
are t, (e 2 ) = A 2 , J', n el = N el /N and n S p = N spm /N. 
N is the size of the system. After normalising by t there 
are two 'coupling constants', A/t and J'/t, and two 'den- 
sities': four dimensionless system parameters in all. We 
will measure all energies in units of t which, from hereon, 
is set to 1. The Coulomb term, even if included, does not 
involve a variable coupling constant. In addition to all 
these there are temperature and magnetic field as exter- 
nal variables. In what follows we will ignore the Coulomb 
term. We will also treat the core spin Si as 'classical', 
and study the magnetic properties, thermodynamics, and 
some aspects of transport in this model. 



III. COMPUTATIONAL SCHEME 



Z = J VS i Tr CtC te-K H " +J "£ i <*- s 4 (3) 

We can formally trace over the fermions and write the 
partition function purely as an integral over spin config- 
urations, i.e, 

Z=fvS i e-^ff^ (4) 
where the effective 'spin Hamiltonian', H e ff, satisfies 

It is obvious that the functional H e f /({Si}) is just the 
fcrmion free energy F e i in the background {Si}. In gen- 
eral it depends on the full spin configuration {Si} and 
not just pairwise interactions. Our problem 'reduces' 
to fermions in the background of some quenched disor- 
der, exchange coupled to spins picked from a distribution 
P({SJ) cx e - pH *tM s A). To calculate physical proper- 
ties involving the fermions, the DOS or conductivity, say, 
we have to average over spin configurations picked from 
the distribution P appropriate to a given temperature. 
The temperature dependence and the spin-spin correla- 
tions implicit in P indicate that the 'disorder' arising 
from the J' term is very different from the quenched un- 
corrected disorder in the e-j's. 

Loosely, the fermion properties can be calculated by 
(weighted) average over probable spin configurations. 
The weights themselves depend on a knowledge of the 
fermion free energy in a given configuration! We have a 
coupled problem here, and need to solve it self consis- 
tently. 



Having defined the Hamiltonian we have to adopt a 
scheme for evaluating its properties. The problem in- 
volves strong disorder, in the variables {ei}, as well as 
temperature (and field) dependent 'disorder' in the spin 
background {Si}. Part of the problem is to determine 
the distribution of spin configurations, appropriate to a 
given T and h, the rest to calculate and average elec- 
tronic properties over the disorder (on site & spin). In 
the regime we are interested in, neither of these tasks is 
analytically tractable. Thankfully, there is an essentially 
exact, and implementable, numerical scheme which can 
be employed here. This is what we discuss next. 

Consider the model for some arbitrary coupling J' and 
disorder A. The spins are considered to be classical (for- 
mally, S — > oo). The problem looks like quadratic (non 
interacting) fermions in the background of a spin config- 
uration {Si}. The task is to determine the appropriate 
background configuration(s), {Si}, at a given tempera- 
ture and specified values of the electronic parameters. 
The problem can be set up formally as follows: Con- 
sider a Hamiltonian H = H e i + J' J2i a i-^i- Here H e i 
is a quadratic fermion Hamiltonian excluding terms in- 
volving the Si's. The partition function of this system 
is: 



A. Exact enumeration 

We have written a formal prescription for the spin dis- 
tribution: how do we proceed any further? A computa- 
tional scheme has been developed and explored ovcpJ-bt: 
last few years for handling this kind of problemsE3'ca 
Starting with some arbitrary configuration {Si}o, the 
Metropolis algorithm is used for updating the spin orien- 
tations and generating new (acceptable) configurations. 
The acceptance or rejection of a spin move, which de- 
pends on nearest neighbour orientations in short range 
spin models, now involves the diagonalisation of a N x N 
fermion Hamiltonian. Unless some approximations can 
be made, simulating such a model involves computa- 
tional effort of 0(N 4 ), where N is the system size. A 
factor of N 3 comes from each diagonalisation (accep- 
tance/rejection of a spin 'move') and another N from 
the need to 'update' all N spins in the system. This 
A 4 process has to be repeated ~ 10 3 — 10 4 times at each 
temperature for equilibriation and averaging over equilib- 
rium configurations. At equilibrium, the spin configura- 
tions generated are a sampling of P({Si}). Fermion aver- 
ages calculated over these configurations represent equi- 
librium average at that temperature. While this scheme 
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is 'exact', and has no sign problems unlike fermion QMC, 
the cost for handling large systems is still prohibitive. On 
a standard DEC Alpha workstation one can handle sys- 
tems ~4x4x4 i.e N ~ 100. Working out the tempera- 
ture dependence of x, say, for fixed electronic parameters, 
involves about 12 hours of CPU time. The largest sys- 
tem studiedE2l is 6 3 for the double exchange model. To 
get a first impression of these complicated systems, the 
thermodynamics and energetics of various phases can be 
reasonably studied for these small sizes. However, where 
transport studies are of interest, and they are the most 
interesting measurable property in these systems, the ac- 
cessible sizes are much too small for making useful state- 
ments (about (Jdc: say). This is because the mean free 
path l m fp in the low resistivity phases is often larger than 
the system size (l m fp ^ L) and no useful statements can 
be made about the conductvity of the thermodynamic 
system. Some advances have been made by (i) making 
the algorithms more efficient that iV 4 , and (ii) using 'ef- 
fective' spin models for generating P({Si}), rather than 
handling the full H e f /({Sj}). Algorithms in the first cat- 
egory, arc . rather complex and still in an experimental 
phaseEilEa, while for approximations in the second cate- 
gory there is often no obvious small parameter. In this 
study we use an approximate scheme, described in the 
next section, after benchmarking it against small system 
simulations based on the exact scheme. 



B. Perturbative expansion in J' 

While there is no obvious approximation for arbitrary 
values of the quenched disorder A and the exchange cou- 
pling J', the experimentally relevant regime in J' (see 
later) allows a controlled approximation. We have ar- 
gued that P({SJ) cx e -^i({S<}), L et USj f or simplicity, 
look at T = 0, in which case P({S 2 }) cx e^^HS*}). If 
the coupling J' were small, in a sense to be quantified 
soon, we can try expanding the energy E e i({Si}), about 
the J' = limit, in powers of J' . This is just quantum 
mechanical second order perturbation theory, estimating 
the change in ground state energy of the electrons due 
to coupling to an arbitrary spin configuration. It is a 
variant of the textbook 'RKKY' argumentBEfl and we 
briefly describe the steps below. 

The reference system (at J' = 0) is spin degenerate, so 
(erf)o = 0, for fi — x,y,z, and there is no 0{J') contri- 
bution to the energy. The 0(J' 2 ) term is 



A£ ({S 4 }) 



rn^O 



N\\2 



(6) 



The 's are many particle wavefunctions (Slater deter- 
minants) including the effect of on site disorder. The 
label N refers to particle number here (not system size) 
and m is the label for the state. Eg and E m refer to 
(unperturbed) ground state and excited state energies of 
the N electron system. Writing out the energy change; 



(^IJ'E^i-Sil^X^IJ'E^i-Sil^) 



E 

J,2 EE 

E -^ySi-S. 



(^| t 7 i .S i |^)(^| < 7 J -.S J |^) 

E a N -E% 



where 



Jij — J 



: E 

m=£0 



yjV\ 



F N 



E N 



(7) 



(8) 



The y, label in a could be x or y or z (not summed 
over). A large part of AE actually comes from the local, 
Sf, term but that is a constant and does not affect the 
'cost' of rotating a spin with respect to its neighbours. 
In defining the effective spin Hamiltonian we look only 
at the i ^ j terms in AE. 

To calculate J^, suppose the single particle wave- 
functions of the reference (disordered) electron problem 
are of the form 7^10) = AfcJ CT |0) and the inverse is 
c- = B l a ^ acr . Since A is a 'rotation matrix' we will have 
B = A^ 1 = A T *, The intermediate state |\E^) is created 
by a particle- hole creation operator, of the form 7^7/3, 
acting on 1^^) so it is indexed by two labels a,(3,a ^ (3. 
Using this notation the exchange constants are specified 
as 



Jij — 



4 ^ 

ajt/3 



{BfB^BfB? + c.c}(^ 



- £/3 



(9) 



The £'s are single particle energies and the n are Fermi 
factors. It can be easily shown that Jy is J /2 Xij(u = 0) 
where Xij 1S the spin susceptibility (response function) 
of the electron system. For electrons in free space this 
is just the Fourier transform of the standard 'polarisabil- 
ity' Xo(q)- The (inverse) energy scale for \o is ~ N{ep), 
i.e we may write xo(q) = N(ep)f (q) where f is a di- 
mensionless function 0(1). Using this, the typical mag- 
nitude of the (nearest neighbour) exchange constant is 
~ J' 2 N(e F ). 

Using similar arguments on a lattice, the exchange 
coupling between nearest neighbour spins (on the lat- 
tice) will be ~ J' 2 /W, since N(e F ) ~ l/W where W is 
the bandwidth. In the disordered electron system such a 
simple identification is no longer possible but the typical 
nearest neighbour coupling will still be ~ J' 2 N(cf)- 

It is apparent that if J' were small there would be 
a hierarchy of energy scales which may help us in an- 
alytically defining the effective spin Hamiltonian. If 
J' < t then for the 'clean' electron system we will have 
t ^> J' ^S> Jij, while for the disordered problem we will 
have A > t 3>r4' Jij- Experimentally the value of 
J' is estimatedE-3 to be < 50 K. Even with a modest 
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t ~ 1000 K, we would have N(e F )J' ~ 0.05 < 1. We 
will use a slightly larger J', i.e J' ft ~ 0.2, and the small- 
ness of N(cf)J' will allow us to work with a quadratic 
spin Hamiltonian. Higher order couplings, between three 
spins etc, will be smaller by this factor. 

The primary scheme, as outlined at the beginning of 
this section, is to average electronic properties by diago- 
nalising the electronic part with spin configurations cho- 
sen from the appropriate distribution. In our approxi- 
mation, the spin distribution is specified as P oc e~P Hap , 
where H sp — J^t^y Jij^i-Sj , and does not have to be 
numerically evaluated. However, we still have to do a 
Monte Carlo (MC), with long range spin-spin interac- 
tions, to generate the equilibrium spin configurations. 
The problem involves the following steps now, (i) spec- 
ify the exchange constants in terms of the spectrum and 
eigenfunctions of the reference 'site disorder' problem, 
(ii) perform classical Monte Carlo on the spin model to 
generate a set of equilibrium configurations, in effect a 
sampling of P({Si}). (Hi) at a given temperature, af- 
ter the spin system equilibriates, diagonalise the electron 
problem with some of the spin configurations, and aver- 
age the computed correlation functions, (iv) since there 
is quenched disorder in the problem, repeat this entire 
process for several realisations of {e,} and average the 
observables. 

If the 'rotation cost' etc can be described by a 
quadratic Hamiltonian then 'polaronic' effects, which im- 
ply that spin rotation energy depends on the presence or 
absence of an electron nearby, are irrelevant. Polaronic 
effects may be relevant, at large J', even in high car- 
rier density systemsH. The smallness of N(ep) J' ensures 
that such effects are not relevant here. 



IV. RESULTS AND DISCUSSION 

We now present our results, starting with a comparison 
of the 'exact' and 'approximate' schemes for our choice 
of J'. We then discuss results on the d.c conductivity, 
in simple situations where the spin distribution does not 
have to be evaluated through simulation. These data, 
nevertheless, highlight the interplay of structural and 
spin disorder, and provide an estimate of the obtainable 
MR. We finally examine simulation results for magnetic 
moments in a (progressively) disordered electronic back- 
ground and the spin glass signatures therein. 



A. Comparison of the exact and perturbative 
schemes 

Second order perturbation theory in J' should be sen- 
sible as long as N(e F )J' <C 1. We have argued that 
J' is small for the systems of interest. There are two 
consequences of this, (i) It is an enormous technical sim- 
plification because it 'deconvolves' the magnetic problem 



from the electronic one. For a given electron density, and 
a realisation of disorder, we only need to compute the ex- 
changes, Jij, once and then work with the resulting H sp . 
Electronic properties can be calculated by diagonalising 
the fermion Hamiltonian with configurations {S^} gener- 
ated by equilibriating H sp . The computational cost re- 
duces from iV 4 , for the exact scheme, to iV 2 . In this ./V 2 , 
one factor of N is to evaluate the rotation cost for a sin- 
gle spin (in our long range model) another TV to update 
all the spins, (ii) Physically it implies that 'polaronic' 
effects, mentioned earlier are not relevant in this param- 
eter regime. The question of a 'dense polaron system' 
does not arise here. 

The single spin rotation cost is the central quantity 
in the simulation. To emphasise that this energy, at 
small J', can be very accurately calculated from H sp 
we have extensively studied this 'rotation cost' under 
various conditions. Fig.l illustrates some of the cases. 
The energy cost of rotation is the difference between the 
fermion energies E{S}} ~ E{S°}, in the exact scheme, 
or H sp {Sj} — H sp {S®}, in the approximate one. {S°} 
and {S}} are two spin configurations connected by a sin- 
gle spin rotation SSj, say. Fig.l illustrates three cases. 
In each of them a reference configuration {S?} is cho- 
sen and E and H sp are computed. A rotation is made, 
changing {S°} to {Sj} and E and H sp are computed on 
this state. We plot E\ — E and ifL — H® p as a function 
of J' for three reference configurations. The reference 
configuration {S°} is random for the circle and square, 
it is ferromagnetic for the triangle. 
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FIG. 1. Energy cost for rotation of a single spin: exact 
-vs- approximate schemes. The energy difference between 
two spin configurations is computed through direct diagonal- 
isation (symbols) and by using H 3p (lines), as a function of 
J'. The configurations are described in the text. Disorder: 
A = 8, system size 14 x 4 x 4. There is a spin at each site. 
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FIG. 2. Exact simulation -vs- Monte Carlo with H ap : 
clean system. J' = 0.2, h = 10~ 5 . The data is averaged 
over 1000 configurations at each T, after equilibriation over 
1000 steps. Total energy (top panel) and magnetisation (lower 
panel). System size 4 3 . There is a spin at each site: N sp — N. 



The symbols are the 'exact' cost, the quadratic curves 
(firm lines) are the approximate cost. It is obvious that 
even in the worst case the 'exact' and 'approximate' re- 
sults match very well upto J' ~ 0.2. 

Since this configuration based testing cannot be ex- 
haustive, we have run the full simulation, using the ex- 
act scheme, and also done MC with H sp , to check if the 
thermodynamics and magnetism come out similar in the 
two approaches. We have used J' = 0.2 and studied a 
4 3 clean system (Fig. 2) as well as a system with strong 
disorder (Fig. 3). We think the similarity is quite con- 
vincing. We have also studied J' = 1 (not shown) where 
the difference is quite significant. 

The principal conclusion from these tests is that we can 
reliably use an effective quadratic Hamiltonian for the 
spins, and there are no polaronic effects in the parameter 
regime we are in. 



B. Results on transport 



To calculate the conductivity and compare with ex- 
perimental data we would have to do a MC on the spin 
problem, having computed the Jy, and use the Kubo 
formula. The conductivity will have to be averaged over 
equilibrium spin configurations. To compute the MR, 
i.e <r(h), the magnetic problem has to be redone in the 
presence of the field, and the cycle repeated. Even after 




0.0001 



Monte Carlo with H. 
net = 0.1, h = 10 



sp • 

size 



FIG. 3. Exact simulation -vs- 
disordered system. A = 8, J' = 0.2. 
4 3 , N sp — N. The data is averaged over 1000 configurations 
at each T, after equilibriation over 1000 steps. Energy (upper 
panel) and magnetisation (lower panel). 

the approximation we have made, this is a time consum- 
ing process. In this paper we do not discuss those results, 
saving it for the discussion specific to GdSi, but highlight 
some simple limiting cases. 

We want to illustrate the interplay of structural and 
spin disorder, and the difference in a between a spin po- 
larised and random system. 

The conductivity is plotted in arbitrary units, i.e with- 
out putting in e 2 jfi and the lattice parameter. 



1. Effect of on site disorder 

First consider the case with only 'on site' disorder. 
This is the traditional Anderson localisation problem. 
We use the Kubo-Greenwood formulaEj to calculate a. 
The data in Fig. 4 displays the variation in conductivity 
(as cf is varied) and the DOS, with increasing disorder. 
<r(e) is the variation in conductivity as the electron den- 
sity, or ep, is varied. It is not the 'optical' conductivity. 

The inset in Fig. 4 shows the variation in the maximum 
conductivity in the band (a ma x say) which occurs at the 
band center here, with increasing disorder. Note the pre- 
cipitious drop (even in this finite size calculation). The 
most recent numerical calculations^ put the Anderson 
transition in the 3d tight binding model at A <~ 16.5. 

To study localisation near the band edge we need to 
follow the size dependence of the conductivity for a fixed 
e (or electron density). If states at a certain energy, e, 
are localised then er(e : L) will decrease with increasing 
L (exponentially at large L). If the states are extended 
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FIG. 4. Effect of on site disorder: conductivity (upper 
panel) and DOS (lower panel) in the presence of site disorder. 
System size 6 x 6 x 16. The inset shows a m ax the conductivity 
at the band center as a function of increasing disorder; notice 
that the a scale is logarithmic. The data is averaged over 
200 — 600 realisations of disorder depending on A. 

then cr(e : L) will tend to a finite asymptote as L — > oo. 
Fig. 5 shows a(e : L) for system size 6 x 6 x L, with L = 
16 and 32, for two strengths of disorder. The crossing 
point is a crude measure of the mobility edge separating 
extended and localised states. cr(e : L) is approximately 




FIG. 5. System size dependence of the 'conductivity', at 
two values of disorder, A = 2 and 8. The system sizes are 
6x6x16 and 6 x 6 x 32. Average over 200 — 600 realisations of 
disorder depending on system size and A. Notice the crossing 
of the curves in each panel. 



independent of L for states in the center of the band, 
implying that our system size is large enough to provide 
an estimate of the conductivity in the thermodynamic 
limit. 

The purpose of presenting these results on the very 
standard site disorder problem is (i) to validate our con- 
ductivity calculation, and (ii) to make an estimate of the 
fraction of states localised as a function of disorder A. 
As a crude estimate, 15% of the states (in the band tails) 
get localised for A ~ 8. This will be useful for fixing 
parameter values in the GdSi problem. 



2. Effect of 'spin disorder': uncorrelated spins 

Now consider how electron scattering off random mag- 
netic moments affect the conductivity. We consider a 
periodic array of spins with random orientations. This 
mimics an (uncorrelated) paramagnetic state. There is 
no on site disorder. The conductivity, computed by di- 
agonalising the electron problem in these random back- 
ground, is shown in Fig. 6 for three values of J'. The data 
is averaged over 100 spin configurations for each J'. Re- 
member that to compute the conductivity of an electron 
system coupled to spins, the spin distribution has to be 
computed, a non trivial task. Here we just assumed an 
uncorrelated random distribution. This is expected to 
provide an estimate of a in the spin glass phase (with 
h = 0). The inset in the top panel shows the reduction 
in the maximum conductivity with increasing J'. 




FIG. 6. Effect of 'spin disorder': conductivity (upper 
panel) and DOS (lower panel) for electrons coupled to ran- 
domly oriented spins. The exchange coupling is J' = 0.1, 0.2 
and 1.0 and there is a spin at each site. There is no on site 
disorder. Inset in the upper panel shows reduction in band 
center conductivity with increasing J'. System size 6 x 6 x 16. 
Data averaged over 100 spin configurations. 
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A similar problem, with J' ft — > oo, has been studied 
earlieiO in the context of the double exchange model. 
System size dependence of the conductivity in our results 
does not indicate localisation in the band tails (upto er- 
rors ~ 1%). In the J' — * oo problems!! it is known that 
less than 1% of the states are localised in the tail of the 
band. Principally, in the small J' limit in which we have 
calculated the conductivity, electron coupling to 'spin dis- 
order' leads to a finite scattering rate ~ N(cf)J 12 , and a 
small resistivity, and not much by way of localisation. 

We suspect that the mean free path at J' = 0.1 and 
0.2 is probably too large, l m / p > L, to reproduce the 
'infinite volume' conductivity. 



3. On site and 'spin disorder': spin at each site 

The problem relevant for us has both on site and spin 
disorder. Fig. 7 shows the interplay of these two effects. 
The top curve in the upper panel corresponds to on site 
disorder only, A = 8. The lower curve (squares) shows 
the additional effect of random spins, coupled to the elec- 
trons with J' = 0.2. If we focus on the band center, 
for instance, the resistivity increases from ~ 1/0.03 to 
~ 1/0.015, i.e, from ~ 33 to ~ 66 (in arb units). This 
difference, call it Ap sp , arises due to spin disorder. What 
would be the resistivity if only this spin disorder were op- 
erative? This can be seen from the lower panel in Fig. 7. 
The 'spin disorder only' resistivity, p° sp is ~ 0.5 -C Ap sp , 
which is two orders of magnitude larger. 



0.03 




FIG. 7. Interplay of on site disorder and spin disorder: 
conductivity in the presence of a random on site potential 
and randomly oriented spins on a lattice (N sp = N). Top 
panel shows data with on site disorder only (A = 8 circles), 
and on site + spin disorder (A = 8, J' = 0.2, squares). The 
bottom panel shows a for J' = 0.2, as in Fig. 6, without any 
on site disorder. 



We think that p Q sp is being underestimated, since 
Imfp ^ L. Nevertheless, even with a factor of 2 error, 
the difference between p° sp and Ap sp is still significant. 

In this disorder regime, the 'on site' and 'spin disor- 
der' effects are far from additive: there is no Mathiessen's 
rule. This would of course be obvious close to the mo- 
bility edge, where spin disorder will lead to localisation 
of additional states, but it is also prominently visible at 
the band center. 

Within the accuracy of our calculation we have not 
been able to see the shift in the mobility edge on intro- 
ducing spin disorder (the shift, again, would be < 1%). 
However, as noted, even in the regime of extended states, 
there is a dramatic increase in resistivity on introducing 
randomly oriented spins, quite out of proportion to the 
small coupling, J' — 0.2. 

The change in a (circles to squares, Fig. 7) is, crudely, 
the difference between YSi and GdSi: i.e the effect of 
on site disorder -vs- on site + spin disorder. This also 
provides an impression of the magnetoconductance ob- 
tainable, at large field. The J' = curve on the top 
panel is also the conductivity for fully polarised moments 
(in which case the coupling becomes irrelevant). 



4- On site and 'spin disorder': dilute random array of spins 

Now for the final degree of realism. Consider the ef- 
fect of randomly locating the spins (as in any amorphous 
structure) together with random orientation, Fig. 8. 
There is also strong structural disorder in the back- 
ground, A = 8. There are three sources of randomness in 
this problem, in the site disorder {e^}, the spin locations 
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FIG. 8. Interplay of on site disorder and spin disorder: 
conductivity in the presence of a random on site potential 
and randomly located array of spins. The fraction of sites 
with spins is n sp = 0.3, the disorder is A = 8 and the system 
size is 6 x 6 x 16. The bottom curve (squares) corresponds 
to random orientation of the spins, one of the top curves (cir- 
cles) corresponds to fully polarised spins, both with J' = 0.2. 
The third set (triangles) corresponds to on site disorder only 
(A = 8, J' = 0) and is displayed for comparison. 
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{Ri}, and the orientations {S^}. The {R^j's were on a 
lattice in Fig. 7. The data reveals that the effect of ran- 
domly locating the moments is of no consequence when 
the moments are polarised, o looks as if there are no 
moments at all (top two curves). This is because of the 
strong structural disorder already present. Polarised mo- 
ments only 'renormalise' the structural disorder from A 
to approximately V A 2 + J 12 m A. Disorder in the orien- 
tation however has a strong effect, just as for periodically 
located spins, and all our arguments about MR, in the 
last section, are applicable here as well. The difference 
between the top curve(s) and the lower curve would be 
the 'magnetoconductance' at strong field. We are doing 
a more extensive calculation to quantify the localisation 
effects near the mobility edge. 



C. Magnetic moments in the electron system: effect 
of increasing disorder 

The previous section set out some of the limits for the 
conductivity, arising from site disorder and randomly ori- 
ented or fully polarised spins. For these results to have 
any relevance the actual magnetic state in the system 
should bear some resemblance to the random magnetic 
state we used. In this section we study the magnetic 
state that arises for moments coupled to electrons in a 
disordered environment, and the conditions under which 
a 'spin glass' state can arise. We first consider a periodic 
system, i. e one with a spin on each site of the lattice and 
study the thermodynamic and magnetic properties as the 
disorder, A, is increased. 
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FIG. 9. Distribution of bonds: nearest neighbour, Pi (J), 
next neighbour, P2(J), and third neighbour, Pz(J), as a func- 
tion of increasing disorder. The spins are on a 6 3 lattice, the 
electron density is 0.1 and J' = 0.2. The spikes (S functions) 
correspond to the clean case (A = 0). The bonds are cal- 
culated for a single realisation of disorder for each A. The 
normalising factor for all the distributions is 4 x 10 . 

The clean system is a model case of 'RKKY', with 
some difference arising from the finite size and 'lattice' 
nature of the system. The strongly disordered syste»s 
also seem to have reasonably long range interactions!^, 
but frustrating instead of regular as in the RKKY case. 

The next figure, Fig. 10, shows the result of simulations 
on H sp . The internal energy E(T) and the weak field 
magnetisation m(T) are monitored. The data is averaged 



1. Periodic array of spins 



The magnetic properties are studied via Monte Carlo 
on the Hamiltonian H sp = £w-y JijSi-Sj as discussed 
earlier. The Jij depend on electron density, J', and the 
specific realisation of disorder in e.;'s. For a clean system, 
A = 0, the bonds Jij depend only on the separation i — j 
irrespective of the location of the two sites. More con- 
cretely, when using periodic boundary conditions (PBC), 
all nearest neighbour bonds will have the same value Ji, 
next nearest neighbour bonds will all be J2, and so on. 

The bond distribution for nearest neighbour bonds, call 
it Pi (J) will be a 8 function. Similarly P2(J), for next 
nearest neighbours, P${J), for third neighbours, will all 
be 5 functions. These bond distributions broaden out 
when the electron system is disordered. 

The bond distribution, upto third neighbour, is shown 
on Fig. 9 for a 6 3 system, with PBC, for A = 2, 4, and 8. 
The 'spikes' are the 'distribution' for A = 0. Notice that 
although the first neighbour coupling remains ferromag- 
netic from A = to 8, the farther neighbour couplings 
are antiferromagnetic and not much weaker than the first 
neighbour coupling (at strong disorder). 
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FIG. 10. Susceptibility, specific heat, and structure fac- 
tor in the bond disordered periodic magnetic system. Upper 
panel, susceptibility (x), lower panel, specific heat (Cv), for 
increasing disorder, A = 0, 2, 4. The system size 6 3 , J' = 0.2, 
and n e i — 0.1. The upper inset shows \ for A = 8, and the 
lower panel shows the evolution of the peak amplitude in the 
structure factor, S(q), with increasing disorder. 
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over ~ 1000 equilibrium configurations at each tempera- 
ture and then over 2 — 8 realisations of disorder, depend- 
ing on A. The following are the important features in the 
data: (i) All the way from A = to A = 8 there is a cusp 
in the susceptibility. The temperature at which it occurs 
reduces from ~ 1.8 x 10~ 4 at A = 0, to 0.5 x 10~ 4 at 
A = 8 (see inset in top panel, Fig. 10). (ii) The magnetic 
specific heat has a well defined peak at A = 0, which gets 
broadened with increasing disorder. The specific heat is 
larger in the bond disordered systems at lower temper- 
ature. It is useful to remember that these are classical 
spins, so the Cy remains finite at T = 0, (Hi) the feature 
in x and Cy suggests an ordering transition in the clean 
system, while the nature of the low temperature state 
in the disordered systems is not obvious purely from Cy 
and x- Since the spins are on a periodic structure, we cal- 
culated the structure factor 5(q) from the ground state 
configuration in each case (inset, lower panel). (Hi) The 
clean system has a single peak in 5(q) at q = 0, 0, 7r, of 
magnitude ~ 0(N 2 ). The system is ordered as expected. 
A reduced peak survives at the same q for A = 2, see 
Fig. 10, but disappears at larger disorder. The feature 
in X? along with the absence of any long range order, 
suggests that the strongly disordered systems are actu- 
ally spin glasses. We are in the process of computing 
the Edwards- Anderson order parameter and checking the 
system size dependence of our results on x and S(q). 



2. Dilute random array of spins 

Finally, let us look at the magnetism in a site diluted 
clean system. This would be the canonical RKKY model 
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FIG. 11. Susceptibility and specific heat in the clean di- 
luted spin system. The spins are distributed on randomly 
chosen 30% of the sites in a 8 3 system. The electron density 
is 0.15, J" = 0.2 and A = 0. Applied field h = 10~ 5 . The 
data is averaged over 1000 MC steps at each temperature after 
equilibriation over 1000 steps. 



except that the real exchange in a lattice model looks 
somewhat different from the simple r~ 3 cos(2kF r ) form 
usually assumed in simulating such models. We show 
the x and Cy in such a system, an 8 3 lattice with spins 
on 30% of the sites, x has the same cusp, as in Fig. 10, 
flattening at low temperature, probably indicative of spin 
glass freezing. We are studying the spin glass order pa- 
rameter for this problem, as well as the effect of on site 
disorder. There has been_satne work on RKKY spin 
glasses in three dimensionEjEj but the results about a 
spin glass transition (in the absence of anisotropy) are 
still inconclusive. There is no work, to the extent we 
know, on models with site dilution as well as disorder. 



V. CONCLUSIONS 

We wanted to explore the interplay of structural disor- 
der and doped magnetic moments in an electron system 
within a simple model. Though we are far from demon- 
strating the experimentally observed features near the 
i-m-t, or detailed properties of the observed spin glass, 
we have highlighted how strong disorder can generate 
most of these effects. For instance, relatively weak 'spin 
disorder' [i.e. small J'), acting on top of structural disor- 
der, can have remarkable effects on transport. This spin 
disorder acting on its own, in a high density electron sys- 
tem, would have led only to weak scattering, masked by 
electron-phonon effects etc. Here the change Ap is large, 
and so is the MR. Similarly, the occurence of a disordered 
magnetic state, which is the key to the i-m-t and MR, is 
again facilitated by structural disorder. We saw how this 
happens, due to frustration, for a periodic array of spins. 
We have not shown results on the actual temperature 
and field dependence of transport, we will discuss these 
separately^. Most importantly, we have not provided a 
hint of how Coulomb effects can be incorporated. Let us 
provide a qualitative discussion. _ 

A controlled perturbative schemecl for handling e-e ef- 
fects exists at weak disorder. This should allow an un- 
derstanding of the metallic phase, far from the transition, 
in terms of the spectrum and wave] 'unctions computed in 
our scheme. Similarly, deep in the insulator, a knowl- 
edge of the localisation length and screening allows us 
to quantify the effects of e-e interactionsEJ. There is no 
detailed theory which allows us to-access physical proper- 
ties near the i-m-t (see the workEJ by V. Dobrosavljevic 
and G. Kotliar though), although there are RG calcu- 
lations (reviewedcl by Belite,and Kirkpatrick). There is 
also an interesting schemec3 suggested by Kuchinskii et 
al. but the 'control' in their procedure is not obvious. 
So, although it seems difficult to start completely from 
first principles, including e-e and electron-spin interac- 
tions, an understanding of our simplified model can help 
in quantifying the e-e effects in the problem. 
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